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Abstract 

The concept of a digital cavity is presented. The functionahty of a 
tunable radio- frequency/microwave cavity with unrestricted Q-factor 
is implemented. The theoretical aspects of the cavity and its potential 
applications in high resolution spectroscopy and synchronization of 
clocks together with examples in signal processing and data acquisition 
are discussed. 

1 Introduction 

The development of analog to digital converters (ADC) is advancing rapidly. 
The cutting-edge ADCs available commercially have a bandwidth of 20 GHz 
and can digitize signals at astonishing speed of about 65 GS/s (giga samples 
per second). This rapid development in ADCs have enabled analysis of the 
signals from DC to microwave frequencies digitally. 

Traditionally resonator circuits (LC circuits) and cavity resonators have 
been used to analyze high frequency signals. The resonance frequency of an 
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LC circuit can be tuned, which makes them useful in apphcations where the 
receiver needs the flexibihty to analyze signals coming at different frequencies. 
However, due to the inherent power losses in the circuit the quality factor 
(Q-factor) that can be achieved in an LC circuit is only about 10^. Radio 
frequency cavities or microwave cavities, made up of a closed metal structure 
that confine electromagnetic field can have extremely high Q-factor (upto 
10^). Each cavity has a set of resonant modes (normal modes) that are 
determined by the physical parameters of the cavity. These are not easy to 
tune without compromising the performance of the cavity. Digital cavities 
overcome the limitation of these resonators, i.e. they are tunable, and at the 
same time, can have extremely high Q-factor. 

In the simple case of a linear cavity formed by two reflecting surfaces, the 
mode spacing or the free spectral range fsr is given by 

/- = ^. (1) 

where c is the speed of light and L is the distance between the surfaces. The 
frequencies of the normal modes are given by N ■ fsr with G {1, 2, 3, ..}. 
A sinusoidal signal at the resonance frequency is amplified by the cavity due 
to constructive interference in the successive round trips in the cavity while 
the other signals are suppressed due to destructive interference. A time- 
periodic boundary condition applied to a sinusoidal signal leads to similar 
interference effects mimicking the functionality of a cavity. Here, we use 
the term time-periodic boundary condition to emphasize that the condition 
applied is similar to the periodic boundary condition (in spatial dimensions) 
used for example in molecular dynamics simulations [H |2l |3] to simulate 
properties of a large system by tessellation of a smaller simulation box. The 
periodic-boundary condition in time is implemented by translating the signal 
x{t) at time t to x{t + k ■ T), where T is the period and k G {0, 1,2, ...}. 
Implementation of such a boundary condition during the digitization of a 
signal by an ADC leads to the formation of a digital cavity. 

Fig{T] shows a pictorial representation of a digital cavity implemented 
using discrete values of digitized signals. In the figure the signal is discretized 
at a fixed interval of At. A period of the digital cavity contains n = 10 data 
points with period T = n ■ At. Sinusoidal signals that are in resonance and 
off-resonance with the cavity are shown in panel (a) and (b), respectively. 
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Figure 1: Digital cavity with n = 10 and cavity length T = 10 At. In (a) the 
signal is resonant with the cavity while it is non-resonant in (b). 

2 Theory 

A digital cavity of length n is defined as 

yij-At;n) := ■ At + /c ■ T);j = l,2,...n, (2) 

k=0 

where the waveform y is the response of the cavity to the signal x and Nc + 1 
is the number of times the cavity function is applied to the signal. The cavity 
defined by Equ.dJ]) amplifies the resonant frequencies while suppressing (or 
filtering out) the non-resonant frequencies. 

The fundamental mode of the cavity is given by /o = l/(n ■ At), which 
can be tuned either by changing n or At. In an analogous way a physical 
cavity can be tuned either by changing the cavity length or the dielectric 
constant of the medium in the cavity. The free spectral range of a digital 
cavity is the same as the fundamental mode {fsr = /o), hence the cavity 
amplifies a frequency comb [1] with frequencies that are integer multiple of 
the fundamental mode, i.e. /o, 2/o, 3/o, ... including the DC signal as shown 
in FigJ21 We denote the resonant frequencies by fres- Note that a DC signal 
is also amplified by an analog cavity. The DC offset of a signal in a digital 
cavity can be filtered out using a high band pass filter or an AC coupling 
prior to digitization. Another simple approach to filter out the DC offset is 
to use a generalized digital cavity as defined later. 

The upper limit of the response of the cavity to a non-resonant frequency 
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Figure 2: A frequency comb that is amplified by a digital cavity. The x-axis 
is presented as frequency coordinate relative to the fundamental mode. The 
amplification as well as the line width of the resonant modes depend on Nc- 



f is given by 



V5 fr 



(3) 



27r / fres 

where fres is the resonant frequency closest to / (see Appendix for the deriva- 
tion), and the norm in the left side of the equation is taken to be the absolute 
value of the waveform. As shown by Equ.Q there is no upper bound for res- 
onant frequencies. In this case the waveform gets amplified linearly with Nc, 
i.e. \\y{j ■ At) 1 1 = {Nc + l)||a;(j ■ At)||, and the cavity response can be made 
arbitrarily large by choosing Nc appropriately. If we define the amplitude 
of the response of the cavity to a signal as A{f) = 'max{\\y{j ■ At)\\) then 
the ratio of the amplitude for a resonant frequency to a non-resonant fre- 
quency, A{fres)/A{f), can be arbitrarily large. This, in other words, means 
that the line-width of the cavity can be arbitrarily small. We take full-width 
at half the maximum (FWHM) as the measure of the line-width. Let C,fo 
denote the half- width at half the maximum. Then for / = /o + ^/o we have 
A{f) = {Nc + l)/2, and using Equ.© we get 
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As the waveform generated after the cavity, y{j ■ At), is a set of numbers, the 
definition of the amplitude can be generalized. If the amplitude is defined 
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Figure 3: Line-width of a digital cavity using r = 2 and different values of 
Nc- The line- width narrows progressively with increasing Nc- The figure 
also shows that at certain non-resonant frequencies, indicated as nodes, the 
cavity completely suppresses the signal. 



as A{f) = maxdKy^j ■ A))*^'!!), where r > 1, then the generalized FWHM is 
given by 

- 7r(iV,+ l)- 

As shown by Equ.([5]) the line-width of a digital cavity can be made arbitrarily 
narrow by choosing large Nc, and in addition it can also be digitally narrowed 
by simply choosing larger r. Figj3] shows line-widths for a digital cavity with 
different values of Nc and r = 2. 

Figj3] also shows nodes at certain frequencies at which a digital cavity 
completely suppresses the signal even for low values of Nc- The frequencies 
at which these nodes occur depend on Nc, and they occur at / = fres + 
Pfo/ i^c + 1) for j9 G {1, 2, Nc}. As pointed out later these nodes can be 
useful in demultiplexing a signal that is multiplexed over different frequencies. 



3 Discussion on the Applications and Scope 

One of the scientific applications of a digital cavity is in phase-sensitive signal 
detection. We developed the concept for real time analysis of phase and am- 
plitude of different frequencies present in the fluorescence signal from a sam- 
ple that is excited with phase modulated laser pulses as in the fluorescence 
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detected wave-packet interferometry [51 E] and the two-dimensional electronic 
coherence spectroscopy. [7] In these experiments, the fluorescence from the 
sample is multiplexed over many frequencies due to the contribution of the 
linear and the non-linear response of the sample. The signals at different 
frequencies provide a plethora of information about the system being investi- 
gated. A data-acquisition system that is capable of de-multiplexing the dif- 
ferent contributions simultaneously would make this technique promising for 
chemical imaging. This in conjunction with near-field scanning techniques, 
can provide detailed information of chemical composition of a system with 
sub-wavelength resolution. [8] A digital cavity can be used to de-multiplex 
the signal and analyze the different contributions simultaneously. Thus, it 
could enable multiplexed chemical imaging using the fluorescence detected 
interferometric techniques. FigJH shows two signals at 4 and 6 kHz (panels 
(b) and (c), respectively) recovered from a signal that is synthesized by dig- 
itally mixing four frequencies at 4, 5, 6 and 7 kHz. The signal is digitized 
at the rate of 96 kS/s (kilo samples per second). Two digital cavities with 
n = 24 and 16 are implemented in real time to recover the two different 
contributions to the signal. Other contributions can also be recovered at 
the same time simply by implementing other cavities with different cavity 
lengths. 

De-multiplexing technique as shown in FigJH can also be used in data 
transmission employing multiple frequencies. When the frequencies are well 
defined, even small values of Nc can select the signal contribution from the 
desired frequency. The 4 kHz signal in FigH] is selected by using Nc = 

3 and the 6 kHz signal is selected by using Nc = 5. Note that for the 

4 kHz signal all the other frequencies, 5,6 7 kHz, lie on the nodes of the 
cavity response when Nc = 3, i.e. they are completely suppressed, and 
the same is true for the 6 kHz contribution when Nc = 5. Multiplexed 
data transmission can dramatically improve the data transfer rate compared 
to the transmission using single frequency as is commonly done in wireless 
communication. Clever implementation of digital cavities can have significant 
impact in this field. 

A digital cavity described by Equ.(|2]) is generalization of coherent sam- 
pling technique [9l [10] in which a common reference clock is used for signal 
generation as well as for digitization. Coherent sampling has been used for 
phase-sensitive detection, similar to that shown in FigJU of signals mainly 
at /o and 2/o to investigate rare gases using high-resolution tunable diode 
laser spectroscopy. [lOj A digital cavity is generalization of the previously 
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Figure 4: Signal de-multiplexing using a digital cavity. Panel (a) shows the 
signal that is generated using direct digital synthesis by summing signals at 
frequencies 4, 5, 6 and 7 kHz. All the frequencies have same amplitude and 
zero phase offset. The signal in (a) is digitized using an ADC at a rate of 
96 kS/s. A digital cavity with n = 24 selects the signal at 4 kHz (panel (b)) 
and another cavity with n = 16 selects the signal at 6 kHz (panel(c)). Nc 
is set to 3 and 5 in (b) and (c), respectively. All the waveforms have been 
normalized by the maximum value. 
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implemented techniques in that it does not need the same reference clock for 
signal generation and digitization. Furthermore the definition of the digital 
cavity can be generalized as 

y{j ■ At) := ^ x{j ■ At + n- At-k) cos{k ■ 6) + 

k=0 

^x{{j + 1) ■ At + n- At-k)sm{k-9), (6) 

k=0 

where / G {n/4, n/2, 1, 2, 3, ..} is used to offset the fundamental frequency /o 
to fg = ra/o/(4/), and 9 = 27rA;(4Z — n)/(4/) (see Appendix for the derivation). 
Equ.(|2]) is a special case of the generalized cavity in which 41 = n . For 
4/ = 2n the generalized cavity functions as a comb-pass filter amplifying 
signals at /o/2, 3/o/2, 5/o/2, and for other / it acts as a comb-pass filter 
amplifying signals at fg,fg + fo, fg + 2/o,.... Note that a DC offset of a 
signal in the generalized cavity, except for the case when 41 = n, becomes 
non-resonant with the cavity. The generalized digital cavity can be tuned 
by changing either of At, n or /. Such tunability significantly increases the 
scope of a digital cavity compared to previously known techniques like the 
coherent sampling. 

As a digital cavity is freely tunable to signal source (s) at a particular 
frequency it can be used in applications like software defined radio in which 
the tuning of the radio to a certain frequency is done by the software. When 
a signal needs to be transmitted over a noisy channel then the cavity can be 
used to suppress the random (white) noise. Note that the cavity averages 
a cycle of the resonant signal times, which improves the signal to noise 
ratio by a factor of \/Wc for the random noise. Fig J5] shows two sinusoidal sig- 
nals retrieved from their corresponding noisy parent signals. Some examples 
where suppression of white noise becomes important is in signal transmission 
over long distance (satellite communication, interplanetary communication, 
etc.), detection of faint radio/microwave frequencies from astronomical bod- 
ies, detection of cosmic microwave background, and perhaps even signals from 
extra-terrestrial intelligent life forms as sought after by the SETI project. In 
some cases deliberate adulteration of signal with white noise can be useful 
to avoid eavesdropping. Digital cavities can significantly improve the instru- 
mentation and reduce the cost related to high-end data acquisition/analysis 
systems needed in the works mentioned above. 
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Figure 5: Reduction of random (white) noise in a digital cavity due to wave- 
form averaging. The signals are generated by mixing random noise five and 
ten times the amplitude of a sine function at 960 Hz in (a) and (b), respec- 
tively. Digitization is done at the rate of 96 kS/s. x— axis shows the index j 
used to denote the digitization position - one cycle is completed for j = 100; 
two cycles are shown in the figures. The signal to noise ratio improves in 
both the cases ((a) and (b)) with larger N^. The dotted lines are fits to the 
cavity response using sine functions. 
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Another potentially important application of the digital cavity can be in 
high resolution spectroscopy. As shown in Equ.([S]) the line- width of the cavity 
becomes narrower with increasing Nc, i.e. data acquisition time. This is the 
consequence of the fact that the precision of any interferometer depends on 
the number of the cycles of the sinusoidal signal that interfere. [11] In analog 
cavities the multi-pass within the cavity increases the number of cycles that 
interfere, however, the power losses in each pass sets an upper limit. Power 
loss in this case equates to information loss. A digital cavity is not prone 
to this limitation as information present in each cycle is stored digitally. 
Consequently, a digital cavity can have arbitrarily high precision, depending 
on the time over which the data is acquired. Data acquired for Is gives 
a precision of about 1 Hz in the measurement of a frequency, and higher 
measurement time results in better precision regardless of the frequency being 
measured. Every digitizer, however, has the problem of electronic jitter, 
which introduces uncertainty in time at which digitizations happen. Jitter 
in digital cavity is analogous to various noise sources in an analog cavity like 
the thermal noise and zero point fluctuation (vacuum fluctuations). In the 
high end digitizers, jitter can be less than 100 fs. The jitter does not set 
an upper limit but increases the data acquisition time to reach the desired 
precision. With the currently available digitizers it would be possible to 
analyze the hyperfine transitions in atoms and molecules in general, and 
atoms used for atomic clocks in particular, with extremely high precision. 
As the digitization is controlled by a clock, precision measurements with a 
digital cavity can also be used to synchronize clocks. When a microwave 
of 20 GHz generated by a source controlled by a master clock is used for 
synchronization then the frequency measurement done by a digital cavity 
can be used to synchronize the clock of the digitizer to ~ 10~^^s, and a day 
of measurement enables synchronization better than 10~^^s , which rivals 
the synchronization that can be achieved using analog techniques based on 
radio/microwave freqeuncies. |1] However, implementing a digital cavity with 
stable operation for a day or longer could be challenging. 



4 Conclusion 

We have presented the theoretical description of the digital cavity and derived 
some of its properties. Digital cavities, in general, implement the functional- 
ity of the analog cavities in time domain and eliminate problems related to 
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power losses in the analog cavities. We have also shown examples of how dig- 
ital cavities can be used in signal analysis and discussed other possible uses in 
signal transmission, high precision spectroscopy and precise synchronization 
of clocks. 
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5 Appendix A: 

Bounded and unboundedness of the cavity response 

To simplify the derivation of the properties of a digital cavity we represent 
the signal x{i ■ At), i G Z as a Fourier series. Without any loss of generality, 
we consider a component of the series, a sine function, sin(a;t -f- digitized 
at time intervals At. We use the following equation to to define a digital 
cavity of length n, n & N"*". 

y{j ■ At; n) := ^ sin(w ■ j ■ At + u ■ k ■ n ■ At + (j)). (7) 

fc=0 

In (I7j) n ■ At = T gives the period of the cavity, thus for each j e (0, 1, . . . , n — 
1) the equation sums up the values at equivalent points in the successive 
periods for Ac times. Instead of using frequency / as done in the main 
text, we use angular frequency u in the equation to reduce notations in the 
equations. We define cUc = 27r/T as the fundamental frequency of the cavity. 
Using this definition, Equ.(I7]) can be written as: 

l/(j-At) = 5^sinfu;-j-At + ^A; + 0j . (8) 

fc=0 ^ J 

We further analyse Equ.® in the following two cases: 
First, when u = Uc, i.e. when the cavity is resonant with the function, 
the cavity simply amplifies the sine function A,, -|- 1 times: 

y{j ■ At) = (Ac + 1) sm{uj, ■ j ■ At + 0). (9) 
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Equ-dH]) is true for any u = mUc , where m G Thus for all the 
frequencies that are integer multiple of the cavity frequency, application of 
the digital cavity amplifies the signal and at the same time increases the 
signal to noise ratio. As there are no any physical constrains that limit the 
amplification, it does not saturate. In other words, the amplification factor 
is not bounded. 

Now we analyze the case when u ^ muc, i.e. when the cavity is not 
resonant with the function. For simplicity we analyse the case Uc ^ u) ^ 
1.5uc- By symmetry, the following argument also applies for the other cases. 
We can rewrite Equ.® as 



A:=0 



. f2nu; 



sm 



k + 



(10) 



where 02 = ■ J ■ At + 0. Let uj/uJc= 1 + x with ^ x ^ 0.5. Equ. ( ITOl) can 
be written as 



y{j ■ At) = cos 02 sin(f2/c) + sin 02 cos{Qk), with Q = 2'kx. 



k=0 



k=0 



The sums on the right hand side have an upper bound, which can be calcu- 
lated as follows. Assume Nc is very large. Take r such that Qr ^ n ^ Q{r+1). 
Then 



sin(n/c) 



fc=0 



k=Q 



(11) 



Let Qk = y, then Ak = Using Ak in the above 



Similarly 



sin(nA;) 



k=0 



E^=o sin y Ay 

nKk 



7r(w - Wc)' 



cos(r2A;) 



k=0 



2n{uj — Uc) 
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Thus for (jj 7^ Uc, J^{j ■ At) is bounded by 

At)|| ^ ^(2cos02 + sin02). (12) 

Note that the singularity at Uc indicates the divergence of digital cavity map 
at Uc, which is also clear in Eq.(j9]). 

The maximum value (2 cos 02 + sin 02) can have is a/5, which gives 

WHj-m^^^^- (13) 

Equ. (fT3ll shows that for the case when u ^ Uc, or in general when u ^ muc, 
the cavity response has an upper bound. 

Determination of frequencies completely suppressed by a digital 
cavity. 

If we take u = Uc/Nc in Equ. ffTOj) then 



k=0 



y{j-At) = $^sin(^ + 02) =0 (14) 



for all j as the summation is over a period of the sine function. Thus for 
a digital cavity in which the successive periods are summed A^^^ + 1 times 
there are Nc frequencies between and Uc given by a; = puc/ (A^c + 1) where 
p e Zandp ^ Nc that are completely filtered out. We define Uc/Nc as the 
node distance. The nodes are present at u = i ■ Uc i: puJc/{Nc + 1) for all 
t e NApG {l,2,...,Arj. 



6 Appendix B: 

Frequencies selected by digital cavity and tuning the cavity. 

According to Equ.([H]) for all u = iuc where i is an integer we have 

y{j ■ At) = {Nc + 1) sin(a; ■ j ■ At + 0). (15) 

Thus the cavity selects all the frequencies that are integer multiple of the 
fundamental frequency, in other words it selects a frequency comb with carrier 
offset frequency and comb tooth spacing both given by Uc- The fundamental 
frequency, and thereby the comb, selected by the cavity can be changed either 
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by changing n (number of samples digitized per unit time) or At (digitization 
interval); remember ujc = ^tt / {n ■ At). However, the definition of the digital 
cavity can be generalized such that carrier offset frequency and comb tooth 
spacing are independent. For e.g. if we define the cavity as 

y{j-At) := ^(-l)'=sin(a;-j-At+^A: + (16) 

then the carrier offset frequency is given by uJc/2 while the tooth spacing is 
given by Uc- Equ ( !T6|l is a special case of the general digital cavity defined 
as 

yU ■ At; 9) := 5^ sin (u ■ j ■ At + 2n (^^^^^^ ^ + j • (17) 
For u to be resonant with the cavity we need 

ZTT -1 

= 2.(i-^) (18) 

The values 9 can have is contrained due to the fact that data is accumulated 
at discrete time points and (f) is not known in general. The general expression 
for the allowed values of 6' in a digital cavity can be found as follows: 
Equ ffTTl) can be written as 



y{i-At) = ^sin (u ■ j ■ At+ '^^k + j cos(6'A;)+^ sin (u ■ j ■ At+ '^^k + + - ) sin(6'A;). 

j,_n V / , _„ \ / 



k=0 ^ ^ / ^=0 

Using the relations 

n ■ At 

and 



71 LOc - n ■ At 



2 

we analyze the second term 



2nd term = sin {uj ■ j ■ At + u ■ n ■ At ■ k + — ■ At ■ Uc + 



k=0 
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Let oOc = ruj, then we have 



2nd term = sin (^u ■ j ■ At + u ■ n ■ At ■ k + — ■ At ■ r ■ u + (pj 



k=0 

According to Equ. ffTSj) we also have 



e = 27v{l--] (19) 



r 



Let 

n ■ r 



4 



/, with / e Z then 



r = ^i. and 

n 

e = 2tt 



4/ 

The frequency that resonates with the cavity is then given by 

- = if- P"' 

Finally the general expression for the cavity can be written as 



y{j -At) = ^ sin(w ■ j ■ At + u ■ n ■ At ■ k + (p) cos (^271 (^—^j—^ + 

^ sin(a; ■ {j + I) ■ At + uj ■ n ■ At ■ k + (p) sin (^27i (—^j—^ ^) 

References 

[1] A. R. Leach, Molecular modelling: Principles and applications, 2nd ed., 
Pearson Education Limited: Harlow, 2001. 

[2] M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids Claren- 
don Press, Oxford University Press, Walton Street, Oxford 0X2 6 DP, 
1991. 



15 



[3] K. Karki et al., Study of structural and dynamic properties of liquid 
phenyltrimethoxysilane, Phys. Chem. Chem. Phys. 13 (2011) 11864. 

[4] S. T. Cundiff, Femtosecond Optical Frequency Comb: Principle, Opera- 
tion, and Applications Springer, Boston 2005. 

[5] N. F. Scherer et al., Fluorescence- detected wave packet interferome- 
try: Time resolved molecular spectroscopy with sequences of femtosecond 
phase-locked pulses, J. Chem. Phys. 95 (1991) 1487. 

[6] P. F. Tekavec et al.. Wave packet interferometry and quantum state recon- 
struction by acousto-optic phase modulation, J. Chem. Phys. 125 (2006) 
194303. 

[7] P. F. Tekavec et al.. Fluorescence-detected two-dimensional electronic co- 
herence spectroscopy by acousto-optic phase modulation, J. Chem. Phys. 
127 (2007) 214307. 

[8] K. Karki et al.. Pump-probe scanning near field optical microscopy: 
Sub-wavelength resolution chemical imaging and ultrafast local dynamics, 
Appl. Phys. Lett. 100 (2012) 153103. 

[9] P. Heinonen et al.. Periodic interference rejection using coherent sampling 
and waveform estimation, IEEE T. Circuits Syst. (1984) 438. 

[10] M. Andersson et al.. Flexible lock-in detection system based on synchro- 
nized computer plug-in boards applied in sensitive gas spectroscopy. Rev. 
Sci. Instrum. 78 (2007) 113107. 

[11] G. W. Bradberry et al.. Measurement of narrow linewidth with a Fabry- 
Perot interferometer of limited resolution Opt. Commun. 20 (1977) 307. 



16 



